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Abstract 

The  reaction  of  ground-state  NH  with  H2  has  been  studied  in  a  high-temperature  photochemistry  (HTP)  re¬ 
actor.  The  NH(X3I7)  radicals  were  generated  by  the  2-photon  193  nm  photolysis  of  NH3,  following  the  decay 
of  the  originally  produced  NH(A3/7)  radicals.  Laser-induced  fluorescence  on  the  NH(A3/7-X3r  0,0)  tran¬ 
sition  at  336  nm  was  used  to  monitor  the  progress  of  the  reaction.  We  obtained  k  (833-1432  K)  =  3.5  x 
10“11exp(-7758  K/T)  cm3  molecule'1  s~lt  with  ±2cr  precision  limits  varying  from  12  to  33%  and  corre¬ 
sponding  accuracy  levels  from  23  to  39%.  This  result  is  in  excellent  agreement  with  that  of  Rohrig  and  Wagner 
[Proc.  Combust.  Inst.  25  (1994)  975]  and  the  data  sets  can  be  combined  to  yield  k  (833-1685  K)  =  4.4  x 
10-11  exp(-8142  K/T).  Starting  with  this  agreement,  it  is  argued  that  their  rate  coefficients  for  NH  +  CO2 
could  not  be  significantly  in  error  [Proc.  Combust.  Inst.  25  (1994)  975].  This,  combined  with  models  of  several 
combustion  systems,  indicates  that  HNO  +  CO  cannot  be  the  products,  contrary  to  their  suggestion  [Proc.  Com¬ 
bust.  Inst.  25  (1994)  975].  Ab  initio  calculations  have  been  performed  which  confirm  this  conclusion  by  showing 
the  barriers  leading  to  these  products  to  be  too  high  compared  to  the  measured  activation  energies.  The  calcula¬ 
tions  indicate  the  likelihood  of  formation  of  adducts,  of  low  stability.  These  then  may  undergo  further  reactions. 
The  NH  +  H2O  reaction  is  briefly  discussed  and  it  is  similarly  argued  that  HNO  -f  H2  cannot  be  the  products,  as 
had  been  previously  suggested. 

©  2006  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

Knowledge  of  the  kinetics  of  reactions  of,  or  lead¬ 
ing  to,  NH  is  important  in  combustion  research.  For 
the  title  reactions,  apparently  the  only  direct  observa¬ 
tions  are  those  of  Rohrig  and  Wagner,  who  studied 
these  reactions  within  the  1100-2000  K  temperature 
range  behind  shock  waves  [1].  NH  in  its  ground  elec¬ 
tronic  state  (X3X)  was  produced  by  thermal  dissoci¬ 
ation  of  HN3.  Its  concentration  was  monitored  by  cw 
laser  absorption.  For  the  reaction 

NH(X3  E)  +  H2(XlE)  -►  NH2(*2B)  +  H(2S), 

AH298  =50±2kJm°rl  O) 

they  obtained  k\  (1156-1685  K)  =  (1.7  ±  0.8)  x 
1 0~ 10  exp(— 10,100  ±  200)  K/T),  for 

^•HNOfX^  +  COfX'r), 

AH^98  =  33  ±  3  kJ  mol- 1  (2) 

they  determined  k2  (1201-1912  K)  =  (1.7  ±  0.8)  x 
10_n  exp(— 7220  ±  200)  K/T),  and  for 

NH(Ar3r)  +  H20(X‘A) 

^-hno(x1a)  +  h2(x-12;), 

AH29s  =  — 8  ±  3  kJ  mol-1  (3) 

they  measured  *3  (1329-1906  K)  =  (3.3  ±  1.6)  x 
10-11  exp(— 6970  ±  200  K/T).  Here  the  rate  coef¬ 
ficients  are  expressed  in  cm3  molecule-1  s-1.  The 
heats  of  formation  are  taken  from  the  JANAF  ta¬ 
bles  [2],  except  for  NH,  NH2,  and  HNO,  which  are 
from  Anderson  [3,4]. 

For  reaction  (2),  Rohrig  and  Wagner  concluded 
that  HNO  +  CO  are  the  products  of  the  CO2  reac¬ 
tion,  as  all  other  product  channels  they  considered 
are  too  endothermic  in  view  of  the  observed  activa¬ 
tion  energy.  For  reaction  (3),  they  identified  several 
other  energetically  accessible  channels,  but  consid¬ 
ered  the  products  shown  as  the  most  likely.  There 
are,  however,  several  difficulties  with  these  conclu¬ 
sions.  When  either  or  both  of  the  reverse  reactions 
(-2)  or  (-3)  with  iiT)  and  k-^(T)  determined 
from  the  above  expression  and  the  thermochemistry 
are  used  in  models  of  several  combustion  systems, 
global  reaction  rates  are  too  fast  by  one  to  four  or¬ 
ders  of  magnitude;  detailed  analysis  revealed  in  all 
cases  that  the  problem  is  intimately  connected  to  the 
inclusion  of  these  reverse  processes.  These  models 
include:  (i)  that  of  Glarborg  et  al.  [5]  for  NO  reduc¬ 
tion  in  rich  NO/CO/H2  mixtures,  in  two  flow  reactors, 
for  the  1200  to  1800  K  temperature  domain,  (ii)  that 
of  Dagaut  et  al.  [6]  for  such  mixtures  in  a  jet-stirred 


reactor  from  800  to  1400  K,  and  (iii)  Anderson’s  cal¬ 
culations  of  the  Diau  et  al.  [7]  results  on  the  thermal 
reduction  of  NO  by  H2  in  mixtures  of  these  com¬ 
pounds  with  CO  in  a  static  reactor  from  900  to  1225  K 
(see  Appendix  A).  Models  for  the  thickness  of  dark 
zones  in  nitrate  ester  propellants  combustion  are  sim¬ 
ilarly  very  adversely  affected  [8].  All  these  groups 
achieved  good  comparisons  of  models  with  experi¬ 
ments  without  reactions  (-2)  and  (-3),  suggesting 
that  the  products  and/or  rate  coefficients  are  incorrect, 
and  thus  decided  to  omit  them  and  their  assumed  re¬ 
verse  processes. 

It  thus  appears  that  either  the  reported  rate  coef¬ 
ficients  are  incorrect  or  that  the  products  of  the  NH 
reactions  with  CO2  and  H2O  are  not  those  indicated 
by  Eqs.  (2)  and  (3).  In  order  to  decide  between  these 
possibilities  we  have  now  studied  one  of  their  reac¬ 
tions  [1],  i.e.,  that  with  H2,  by  a  very  different  experi¬ 
mental  technique,  and  have  made  an  ab  initio  study  of 
reaction  (2).  The  results  indicate  that  HNO  -I-  CO  and 
HNO  +  H2  are  not  the  products  of  reactions  (2)  and 
(3),  respectively.  Much  further  work  will  be  necessary 
to  identify  the  actual  product  channels. 


2.  The  NH  +  H2  reaction 

2.1.  Experimental  technique 

A  schematic  of  the  HTP  reactor  used  is  shown  in 
Fig.  1.  It  and  the  data  acquisition  procedures  have 
been  described  frequently  [9,10].  Briefly,  the  5  cm 
i.d.,  30  cm  long  ceramic  reaction  tube  is  surrounded 
by  resistance  heating  elements  and  insulation  in  a 
stainless-steel  vacuum  chamber.  The  bottom  (up¬ 
stream)  plate  contains  an  inlet  for  N2  bath  gas  and  a 
water-cooled  inlet  through  which  the  photolyte  NH3 
and  the  second  reactant  gas,  here  H2,  flow  sepa¬ 
rately.  Teledyne-Hastings  mass-flow  controllers  are 
used,  and  pressure  is  measured  with  an  MKS  Bara- 
tron  transducer.  The  linear  gas  velocities,  v,  are  fast 
enough  to  provide  each  photolysis  pulse  with  a  fresh 
reaction  mixture.  Residence  times  between  the  cooled 
inlet  and  the  observed  reaction  zone  are  long  com¬ 
pared  to  the  reaction  times  and  are  chosen  such  that 
mixing  is  at  least  99%  complete  [11,12]. 

The  reactor  contains  four  side  ports  at  right  an¬ 
gles.  Two  facing  ports  have  Brewster  angle  win¬ 
dows  through  one  of  which  the  193  nm  photolysis 
laser  radiation  (Lambda  Physik  Compex  201)  is  in¬ 
troduced.  The  other  window  provides  the  entrance 
port  for  the  collinear  diagnostic  beam  from  a  Lambda 
Physik  LPX  100  excimer  FL  3002  dye  laser  combi¬ 
nation.  The  photolysis  radiation  can  be  focused  on 
the  center  of  the  reaction  tube  with  a  plano-convex 
lens.  The  dye  laser  is  tuned  to  336  nm  to  pump  the 
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Fig.  1 .  Schematic  of  the  HTP  reactor. 
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Summary  of  the  rate  coefficient  measurements  of  the  NHtX3^)  ±  H2  -►  NH2  +  H  reaction 


T 

<K) 

P 

(mbar) 

[M] 

(1018cm"3) 

[HjJmax 

(1016cm"3) 

LP 

(ml) 

ENH31 

<1014  cm"3) 

t 

(cm) 

V 

(cms~l) 

(*i  ±&ki) 

<10— 14  cm3  molecule”1  s”1) 

833 

138 

1.20 

9.71 

75 

5.37 

12 

23. i 

0.47  ±0.11 

869 

219 

I.&2 

6.35 

81 

9.43 

16 

14.9 

0.30  ±0.23 

934 

305 

2.36 

15.40 

111 

6.90 

12  - 

11,5 

0.69  ±0.14 

978 

181 

1,34 

12.20 

24 

5.05 

12 

2L2 

1,07  ±0,17 

1025 

138 

0.97 

8.08 

66 

4.25 

12 

28.4 

2. 16  ±0.35 

1045 

162 

1.13 

7.40 

93 

5.00 

3 

24.2 

2,68  ±0.21 

1113 

313 

2,04 

13.10 

i0i 

5.96 

12 

13.3 

3,38  ±0.53 

1155 

163 

1.02 

7.47 

102 

3.00 

8 

26.4 

3.55  ±0.92 

1174 

144 

0.89 

2.65 

90 

2,73 

16 

30.2 

2.96  ±0  33 

1265 

173 

0.82 

5,28 

60 

3.30 

12 

33.3 

10.60  ±±.20 

1368 

278 

1.47 

6,93 

108 

6.67 

3 

1ST 

11. 90  ±1.30 

1432 

72 

0.36 

2.11 

123 

1.57 

12 

74.2 

16.60  ±2.80 

NH(A3f7-X52?  0,0)  transition.  The  resulting  fluo¬ 
rescence  is  observed  with  an  EMI  9813  QA  photo¬ 
multiplier  tube.  The  PMT  signal  is  processed  with 
Stanford  Research  SR  250  boxcar  integrators,  which 
provide  a  range  of  delay  times  from  the  photolysis 
pulse. 

The  experiments  were  carried  out  under  pseudo- 
first-order  conditions,  with  [NH]  [H2]  <C  [Nj]. 
Under  these  conditions,  the  fluorescence  intensity  /, 
which  is  proportional  to  [NH],  is  expressed  by 

/  =  h  expt-fcpsi  0  +  (4) 

Here  kp$\  is  the  pseudo- first-order  rate  coefficient, 
/q  +  B  is  the  intensity  at  time  t  =  0,  and  B  is  the  back¬ 
ground,  due  mainly  to  scattered  light.  The  values  of 
Jtp3i  were  obtained  by  a  weighted  fit  of  the  observed 
/  vs  t  profiles  to  Eq.  (4)  [13].  Typically,  five  k^s\  mea¬ 
surements  at  varying  [H2]  were  used  to  obtain  the 


bimolecular  rate  coefficients,  ki,  at  the  temperature 
and  pressure  of  the  experiment.  The  minimum  [H2J 
values  were  set  ca.  5  times  lower  than  the  maximum 
values  listed  in  Table  L  The  exponential ity  of  the  / 
vs  t  plots  was  tested  by  a  two-stage  residual  analy¬ 
sts  [14].  The  analysis  consists  of  a  visual  inspection 
of  the  residual  plot  followed  by  the  runs  test.  A  new 
objective  validation  method,  based  on  the  combina¬ 
tion  of  the  x2  goodness-of-fit  and  bootstrap  (similar 
to  the  Monte  Carlo)  method,  is  also  used  [15],  The 
two  methods  have  been  found  to  be  nearly  always  in 
agreement. 

The  gases  used  were  N2  (99.995%  from  the  liquid) 
and  NH3  (99.995%)  from  Praxair  and  H2  (99.999%) 
from  Matheson.  The  H2  was  dried  by  passing  through 
drierite.  The  NH3  was  used  as  0.3%  NH3  in  N2,  pre¬ 
mixed  in  the  laboratory. 
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2.2.  NH  production 

There  are  a  variety  of  means  of  producing  NH 
in  the  literature,  most  of  which  are  not  suitable  for 
photolysis  experiments  at  or  near  combustion  tem¬ 
peratures.  The  method  developed  here  is  based  on 
the  observation,  by  Kenner  et  al.  [16],  that  193  nm 
two-photon  photolysis  of  NH3  proceeds  in  part  as 
NH3  +  hv  NH£  +  H  followed  by  NHj  +  hv 
NH(A3/7)  4-  H.  The  NH(A3/7)  has  a  radiative  life¬ 
time  of  4.2  x  10-?  s  [1 73.  Hence,  on  a  time  scale  short 
compared  to  the  >  10-4  s  of  the  present  experiments, 
the  photolysis  leads  to  NH(X3£).  There  is  also  some 
direct  formation  of  NH(X^),  as  well  as  of  the  sin¬ 
glet  states  NH(tfJA)  and  of  NHffc1  E)  [18].  The  a 
state  has  been  shown  to  be  quenched  by  N2  [19], 
kQ  (300  K)  =  (6-9)  x  t0“14  cm3  molecule”1  s”1, 
which  was  therefore  used  as  the  bath  gas  and  as  the 
NH3  diluent.  The  highest  NH  fluorescence  intensi¬ 
ties  were  obtained  by  slightly  defocusing  the  photol¬ 
ysis  laser  beam.  Apparently  this  reduced  photolysis 
of  NH. 

2.3.  Results  and  discussion 

The  measured  rate  coefficients,  and  the  conditions 
under  which  they  were  obtained,  are  given  in  Ta¬ 
ble  L  To  verify  that  the  rate  coefficients  are  inde¬ 
pendent  of  the  reaction  parameters,  plots  of  residu¬ 
als  [AfT)  -  kf]/k(T)  versus  these  parameters  were 
prepared.  Here  k(T)  represents  the  rate  coefficients 
obtained  from  the  fit  expression  Eq.  (5)  below,  and 
Jfcj-  are  the  individual  rate  coefficients  measured.  The 
parameters  thus  checked  are  as  follows:  P,  total  pres¬ 
sure;  [M],  the  corresponding  total  gas  concentration; 
the  second  reactant  and  photolyte  concentration;  LP> 
the  photolysis  laser  power;  zt  the  cooled  inlet  to  ob¬ 
servation  plane  distance;  v,  the  average  gas  veloc¬ 
ity;  the  corresponding  residence  time  t ;  the  products 
LP[H2]  and  LP{NH3],  and  71-1 .  These  residual  plots 
are  found  to  be  independent  of  these  parameters. 

Fig.  2  shows  the  present  data.  These  may  be  fit¬ 
ted  by  the  Marquardt  algorithm  [20]  to  the  form 
Aexp(“E^/T),  where  and  op/T  —  ±2%  con¬ 
tribute  to  the  weighting  of  each  point.  The  fitted  ex¬ 
pression  is 

*1  (833-1432  K)  =  3.49  x  10-11  exp{— 7758  K/T) 
cm3  molecule-1  s  .  (5) 

The  variances  and  covariances  are  a2  =  1.8332  x 
10”1  A2,  a\  =  2.3096  x  105,  and  <t2£  =  203.6  A. 
The  resulting  2 o  precision  levels  of  the  fit  are  be¬ 
tween  ±12  and  ±33%.  Allowing  liberally  for  ±20% 
systematic  errors  yields  accuracy  levels  from  ±23  to 
±39%. 


T  K 

2QM  1500  ’  1000  000 


Fig.  2,  Arrhenius  plot  of  the  NH  -I-  rate  coefficients.  (O) 
k,  present  work;  (A)  k ,  Rohrig  and  Wagner  [1];  ( — )  best  fit 
to  all  data,  Eq.  (6). 


Also  shown  in  Fig.  2  are  the  1156-1685  K  data 
from  Rohrig  and  Wagner  for  this  reaction.  The  agree¬ 
ment  is  very  good  so  that  the  two  data  sets  can  be 
combined.  This  yields  . 


ki  (833-1685  K)  =  4.34  x  10-n exp(-8142  K/T ) 
cm3  molecule” 1  s” 1 .  (6) 


Those  authors  also  noted  fair  agreement  with  and 
data  obtained  by  fitting  observations  in  more 
complex  environments,  i.e,,  for  reaction  (1)  in  the 
2600-2800  K  range  [21]  and  from  the  reverse  reac¬ 
tion  [22,23]  also  at  higher  temperatures. 


3.  The  NH  +  CO2  reaction 


'  -  j": 

J.  /.  The  rate  coefficients 

. .  * '  .  -4*-  *  * 


The  fact  that  two  very  different  techniques  (HTP 
and  shock  tube)  using  different  NH  production  meth¬ 
ods  (NH3  photolysis  and  HN3  thermal  dissociation) 
yield  such  excellent  agreement  for  it]  (T)  makes  it  un¬ 
likely  that  there  can  be  a  basic  error  in  the  technique 
used  for  the  Rohrig  and  Wagner  measurements  on  re¬ 
action  (2).  It  should  be  considered  if  other  processes 
could  have  significantly  influenced  their  NH*  con¬ 
sumption  measurements.  Conceivably,  O  atoms  from 
CO2  dissociation  could  have  interfered  due  to  the 
very  fast  reaction  of  O  with  NH,  leading  to  NO  + 
H  and  OH  ±  N,  for  which  the  rate  coefficients  are 
near  the  gas  kinetic  collision  rate  coefficients,  i.e,, 
on  the  order  of  1  x  lG^i0  cm3  molecule-1  s-1,  es¬ 
sentially  independent  of  temperature  in  the  range  of 
interest  [24],  However,  extrapolated  CO2  dissociation 
data  [25]  down  to  1912  K  indicate  the  [O]  to  be  too 
small  by  several  orders  of  magnitude  to  have  influ¬ 
enced  their  results  [1].  For  an  impurity  to  have  caused 


A.  Fontijn  er  at,  /  Combustion  and  Flame  145  (2006)  543-551 


541 


the  measured  rates,  these  would  have  to  have  been  in¬ 
troduced  with  the  C02,  as  the  comparison  above  for 
the  H2  reaction  suggests  no  such  effect  on  those  ex¬ 
periments*  The  purity  of  their  C02  was  99.995%  [}]. 
Hence  any  impurity  from  that  source  could  only  have 
led  to  reactions  at  rates  orders  of  magnitude  smaller 
than  those  observed. 

We  also  examined  Rohrig  and  Wagner’s  data 
analysis  using  the  detailed  mechanism  and  technique, 
described  briefly  in  Appendix  A.  To  this  mechanism 
we  added  reaction  (2)  and  fc2(T)  from  their  work,  and 
reactions  for  the  HN3  pyrolysis  chemistry*  namely 


CO  is  the  most  likely  cause  of  the  discrepancy  be¬ 
tween  the  results  of  Ref*  [1]  and  the  combustion  mod¬ 
eling  studies  noted  in  the  Introduction*  To  investigate 
this  further,  ab  initio  studies  have  been  performed 
(see  below).  It  may  also  be  noted  that  we  attempted 
to  measure  fc2(T)  by  the  same  technique  as  used  for 
k\(T),  However,  in  this  case  O  atoms  were  found  to 
be  produced  in  the  photolysis  process*  These  could 
interfere*  and  therefore  this  approach  was  not  pursued 
further. 

3.2.  Ah  initio  calculations  and  discussion 


HN3  +  M-»NH-|-N2+M, 

k  =  1.2  x  10-9exp(— 18,042  K /T)  [26],  (7) 

HN3  +  H-»-NH2  +  N2> 

k  =  3  x  10-11  exp(— 2315  K/T)  [27],  (8) 

HN3  +  OH -►  HNOH  +  N2,  t  =3  x  10-11,  (9) 

HN3  +  NH-+NH2  +  N3,  Jt=2xl0-11,  (10) 

HN3  +  NH->-N2H2  +  N2,  k=2x  10-11.  (II) 

All  these  ks  are  in  cm3  molecule"1  s1,  The  values 
of  Eqs*  (9H11)  are  estimates*  Within  their  quoted 
error  limits  [1],  no  influence  on  fc2{T)  was  found 
from  these  other  reactions.  In  addition,  we  considered 
whether  H20  present  as  contaminant  in  the  shock 
tube  could  have  affected  results*  This  appears  to  be 
the  most  likely  possible  contaminant.  For  this  pur¬ 
pose,  1000  ppm  H20  were  added  to  the  modeled  mix¬ 
tures.  Note  that  this  is  probably  at  least  a  factor  of  100 
larger  than  one  would  likely  encounter  in  such  shock- 
tube  experiments  [28]  *  The  effect  of  the  added  H20  is 
negligible. 

Thus,  since  the  Jt2<r)  data  appear  substantiated, 
the  assignment  of  the  reaction  products  as  HNO  + 


Reactants,  products,  transition  states,  and  po¬ 
tential  intermediates  were  characterized  with  the 
GAUSSIAN  98  program  suite  [29].  First,  geome¬ 
tries  and  vibrational  frequencies  were  derived  at  the 
E3LYP  level  of  hybrid  density  functional  theory,  us¬ 
ing  the  6-3 1  IG(d,  p)  atomic  basis  set.  Second,  single- 
point  energies  were  obtained  by  approximate  extrapo¬ 
lation  of  coupled-cluster  theory  to  the  complete  basis 
set  limit,  via  application  of  the  CBS-QB3  approach 
of  Petersson  and  co-workers  [30].  These  energies,  to¬ 
gether  with  zero-point  vibrational  and  thermal  correc¬ 
tions,  yield  enthalpies  relative  to  the  reactants,  3NH 
-I-  C02,  at  298  K. 

Some  of  the  ab  initio  results  are  summarized  in 
Figs.  3  and  4;  complete  structural  details  are  given  in 
the  Supplementary  material.  Table  2  lists  the  CBS- 
QB3  enthalpies  at  298  K  and  the  enthalpies  rela¬ 
tive  to  reactants  for  all  adduct  species  found;  these 
adducts  are  shown  in  Fig.  4.  The  first  7  structures  are 
singlet  molecules;  ClS/6-311G(d,  p)  calculations  in¬ 
dicate  that  in  each  case  the  lowest  triplet  energy  is 
higher.  Structure  8  of  Fig.  4  is  a  triplet  adduct  between 
NH  and  C02. 


Fig.  3.  Potential  energy  diagram  for  NH  +  C02  showing  paths  to  HNO  +  CO.  Circled  numbers  1  and  7  correspond  to  those 
numbers  In  Table  2  and  Fig.  4. 
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Fig.  4.  Various  HNCO2  adduct  configurations. 


Table  2 

Ab  initio  enthalpies  for  stable  HNCO2  species  derived  via 
the  CBS-QB3  method 


Molecule 

CBS-QB3  enthalpy 
at  298  K  (au)a 

Enthalpy  relative 
to  3NH  +  CO2 
at  298  K  (kJ  mol-1) 

1 

-243.51270 

-8 

2 

-243.48645 

61 

3 

-243.49393 

41 

4 

-243.52232 

-33 

5 

-243.47184 

99 

6 

-243.36024 

392 

7 

-243.46509 

117 

8 

-243.42968 

210 

a  1  au  »  2625.5  kJ  mol  1 . 


Fig.  3  shows  potential  energy  diagrams  for  the 
reaction  of  triplet  and  singlet  NH  with  COj.  Spin- 
allowed  abstraction  to  form  triplet  HNO  is  possible, 
but  there  is  a  large  barrier  of  226  kJ  mol"1.  If  inter- 
system  crossing  to  a  singlet  surface  occurs  early  in  the 
reaction,  then  a  cyclic  adduct,  with  N  attached  across 
a  C-O  bond,  can  form,  shown  as  l  in  Figs.  3  and  4. 
This  adduct  is  too  weakly  bound,  by  only  8  kJ  mol"1 
relative  to  reactants,  to  be  a  final  sink  for  3NH,  but  it 
can  undergo  further  chemistry.  The  NH  group  can  in¬ 
sert  into  the  C-O  bond  over  a  barrier  of  126  kJ  mol" 1 
to  form  another  adduct,  shown  as  7  in  Figs.  3  and  4.  In 
turn,  this  second  adduct  can  rearrange  while  travers¬ 
ing  a  small  barrier  and  dissociate  to  *HNO  4-  CO. 


Both  these  paths  to  HNO  4-  CO  exhibit  major  bar¬ 
riers  in  excess  of  the  endothermicity,  and  therefore 
are  inconsistent  with  the  measured  activation  energy, 
60  kJ  mol"1 ,  corresponding  to  Eq.  (2). 

The  calculations  thus  confirm  the  conclusion  from 
the  combined  experimental  and  modeling  work  that 
HNO  +  CO  cannot  be  the  products  of  the  NH  4-  CO2 
reaction. 

An  alternative  product  path  is  creation  of  a  sta¬ 
ble  or  intermediate  adduct.  We  have  tried  a  variety 
of  initial  geometries  and  from  these  derived  opti¬ 
mized  geometries  of  8  candidates,  which  are  drawn 
in  Fig.  4.  Their  enthalpies  relative  to  3NH  4-  CO2 
are  summarized  in  Table  2.  The  most  stable  is  the 
HONCO  structure  4,  but  even  this  is  only  exothermic 
by  33  kJ  mol"1.  It  is  therefore  insufficiently  stable  to 
act  as  a  permanent  sink  for  NH  at  elevated  tempera¬ 
tures. 

The  known  thermochemistry  in  this  system  pro¬ 
vides  a  check  on  the  accuracy  of  the  ab  initio  calcula¬ 
tions.  The  overall  enthalpy  change  at  298  K  for  HNO 
-f  CO  formation  is  computed  to  be  30  kJ  mol"1,  in 
accord  with  the  thermochemical  value  of  33  kJ  mol"1 . 
As  a  second  check,  the  enthalpy  change  for  3NH  + 
CO2  -►  HCN  4-  O2  is  computed  to  be  166  kJ  mol"1 , 
again  in  good  agreement  with  the  literature-based 
value  [2,3]  of  172  kJmol"1.  Computed  AH^g  val¬ 
ues  therefore  seem  to  be  usefully  reliable.  An  em¬ 
pirical  estimate  of  the  singlet-triplet  gap  in  NH  is 
150  kJmol"1  and  prior  theoretical  work  suggests 
higher  values  of  157  and  170  kJmol"1  [31].  De¬ 
generate  orbitals  make  inappropriate  for  the 
single-reference  treatments  applied  here;  neverthe¬ 
less  the  computed  adiabatic  singlet-triplet  gap  of 
170  kJ  mol"1  is  close  to  these  prior  estimates.  In  the 
case  of  HNO,  the  singlet-triplet  gap  is  predicted  to 
be  77  kJ  mol"1,  which  is  equal  to  the  value  derived 
from  a  chemiluminescence  spectrum  [32].  Somewhat 
lower  theoretical  values  led  Gurvich  et  al.  to  rec¬ 
ommend  72  ±  6  kJmoF1  [33].  The  magnitude  of 
this  gap  makes  the  spin-allowed  path  for  the  NH  4- 
CO2  reaction  going  to  triplet  HNO  4*  CO  over  100 
kJmol-1  endothermic,  which  is  incompatible  with 
the  observed  activation  energies. 


4.  The  NH  4-  HjO  reaction 

While  we  have  not  explicitly  studied  this  reaction 
a  few  remarks  about  it  can  be  made.  Regarding  the 
rate  coefficients,  one  possible  concern  is  that  H2O 
might  have  dissociated  thermally  in  the  shock  experi¬ 
ments.  However,  calculations  based  on  recommended 
rate  coefficients  [34,35]  indicate  the  H2O  dissociation 
products,  H  and  OH,  would  have  been  present  in  far 
too  small  concentrations  to  have  substantially  reduced 
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[NH].  Highly  purified  H2O  was  used  [1]  and  a  major 
impurity  effect  on  the  rate  coefficients  is  again  highly 
unlikely,  A  check  of  the  data  analysis  of  Ref,  [1], 
similar  to  that  performed  for  the  NH  +  CO2  reac¬ 
tion  above,  also  showed  no  significant  influence  on 
*3(7)  by  other  reactions.  Thus,  the  discrepancy  with 
the  models  can  reasonably  be  attributed  to  the  prod¬ 
ucts  not  being  HNO  4-  H2.  In  this  case,  alternate 
energetically  accessible  addition  products,  H3NO  and 
H2NOH,  have  been  suggested  [I],  but  these  possibil¬ 
ities  need  further  study. 


5.  Conclusions 

Using  an  independent  technique,  the  NH  +  Hj 
rate  coefficients  of  Rating  and  Wagner  have  been 
confirmed  and  the  temperature  range  of  the  observa¬ 
tions  has  been  extended.  This,  and  further  arguments, 
leads  to  the  conclusion  that  their  NH  +  CO2  and  NH 
+  HjO  rate  coefficients  are  also  beyond  reasonable 
doubt.  However,  as  discussed  in  the  Introduction  and 
further  in  Appendix  A,  modeling  of  three  experiments 
(Refs.  [5,6]  and  the  present  modeling  of  Ref,  [7])  and 
of  propellant  dark  zones  is  adversely  affected  if  ei¬ 
ther  or  both  of  reactions  (2)  and  (3)  and  their  reverse 
processes  are  included  in  the  mechanisms.  Detailed 
analysis  indicates  in  all  cases  that  this  is  due  to  the 
reverse  reactions  providing  a  radical  source  which 
causes  a  vast  overprediction  of  the  global  chemical 
rates.  Excellent  results  are  achieved  when  these  re¬ 
actions  are  omitted.  These  observations  thus  indicate 
that  the  products  cannot  be  HNO  +  CO  or  HNO  + 
Hj,  respectively. 

Furthermore,  ab  initio  calculations  of  the  poten¬ 
tial  energy  surfaces  for  the  NH  +  CO2  reaction  were 
performed  in  the  present  work.  These  show  that  the 
barriers  leading  to  HNO  -f  CO  products  are  much  too 
high  to  agree  with  the  activation  energy  measured  in 
Ref.  (1].  The  calculations  also  indicate  that  adducts 
can  form.  These  may  undergo  further  reaction.  How¬ 
ever,  since  no  unique  products  have  been  identified, 
other  types  of  experiments  and/or  more  extensive  ab 
initio  calculations  are  desirable. 
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Appendix  A*  Significance  of  reactions  (-2)  and 
(—3) 

As  noted  in  the  Introduction,  inclusion  of  k^2 
and/or  k- 3  can  cause  models  to  strongly  deviate  from 
experiments.  This  may  be  illustrated  by  the  effect  of 
using  k^2  or£_3  in  the  model  for  CO2  production  in 
NO/H2/CO  mixtures,  diluted  in  Ar,  as  measured  by 
Diau  et  al.  [7],  The  base  mechanism  used  was  devel¬ 
oped  for  propellant  combustion  [8],  It  is  based  on  that 
of  Miller  and  Bowman  [36]  for  nitrogen  chemistry 
in  combustion,  as  further  extended  in  Refs,  [37,38], 
The  modeling  utilizes  the  SEN  KIN  time-dependent, 
homogeneous  gas-phase  reaction  code  [39],  The  ther¬ 
mochemistry  is  taken  from  Refs,  [2-4],  Fig,  5a  shows 
that  there  is  good  agreement  between  the  calculated 
[CO2]  using  this  base  model  and  the  experimental  re¬ 
sults,  Conditions  are  0.960  bar,  1000  K,  and  initial 
mole  fractions  of  H2,  NO,  and  CO  of  0.0263, 0.0161, 
and  0,0430,  respectively. 

Fig.  5b  shows  the  effect  of  adding  HNO  +  CO 
NH  +  CO2  to  the  model,  assuming  its  rate  to  be  deter¬ 
mined  via  the  k2  of  Ref.  [1  ]  and  the  thermochemistry. 
The  CO2  production  rate  increases  by  a  factor  104. 
Fig.  5c  shows  the  effect  of  instead  adding  HNO  4-  H2 

NH  4-  H2O  to  the  model,  again  assuming  the  rate 
to  be  determined  via  the  k$  of  Ref.  [1]  and  the  ther¬ 
mochemistry,  The  CO2  production  rate  increases  by 
a  factor  103,  The  error  becomes  cumulatively  worse 
if  the  two  reactions  are  added  together  (not  shown). 
Similar  erroneous  results  are  obtained  for  other  exper¬ 
imental  conditions  of  Diau  et  al.  Therefore,  the  actual 
rates  of  these  two  HNO  reactions  appear  to  be  vanish¬ 
ingly  small. 

The  possibility  that  an  error  in  the  thermodynam¬ 
ics  used  in  reversing  reactions  (—2)  or  (—3)  might 
have  contributed  to  the  overprediction  of  the  global 
rate  was  checked  by  assuming  an  error  of  20  kJ  mol" 1 
in  their  heats  of  reaction.  This  is  at  least  a  factor 
of  5  larger  than  the  error  limits  in  the  heats  of  for¬ 
mation  would  permit.  For  the  conditions  of  Fig,  5, 
the  error  in  either  of  the  computed  reverse  rate  co¬ 
efficients  would  be  about  a  factor  of  10.  Using  this 
factor  in  the  calculation  [inputting  reactions  (2)  and 
(—2)  or  (3)  and  (—3)  as  separate,  one-way  reactions 
to  CHEMKIN  and  reducing  the  rate  coefficient  for 
(-2)  or  (—3)],  the  time  for  CO2  growth  is  increased 
(and,  thus,  the  reaction  rate  decreased)  in  both  cases 
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Fig,  5.  Influence  of  the  assumed  reactions  (—2)  and  (—3). 
Note  the  vastly  differing  time  scales  of  (a)-{c),  (a)  (•)  Ex¬ 
periments  from  Diau  et  aJ.  [7];  ( — )  prediction  from  the  base 
model  (b)  ( — )  Reaction  (—2)  added  to  the  model,  using  Jtj 
from  Rohrig  and  Wagner  [I],  (c)  (— )  Reaction  (—3)  added 
to  the  model,  using  A3  from  Rohrig  and  Wagner  [1], 

by  about  a  factor  of  10.  The  rate  is  still  far  overpre- 
dtcted. 


Supplementary  material 

Supplementary  data  for  this  article  may  be  found 
in  the  online  version  at  DOI:10.1016/],combustflame, 
2005.12,012.  Thble  15  lists  Cartesian  coordinates  for 
the  adducts  and  transition  states  shown  in  Figs.  3 


and  4,  derived  at  the  B3LY P/6-3 1  IG(d,  p)  level  of  the¬ 
ory. 
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